Melting of stripe phases and its signature in the single-particle spectral function 
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Motivated by the recent experimental data [Phys. Rev. B 79, 100502 (2009)] indicat- 
ing the existence of a pure stripe charge order over unprecedently wide temperature range in 
Lai.8-a;Euo.2Sra;Cu04, we investigate the temperature-induced mehing of the metaUic stripe phase. 
In spite of taking into account local dynamic correlations within a real-space dynamical mean-field 
theory of the Hubbard model, we observe a mean-field-like melting of the stripe order irrespective of 
the choice of the next-nearest neighbor hopping. The temperature dependence of the single-particle 
spectral function shows the stripe induced formation of a flat band around the antinodal points 
accompanied by the opening a gap in the nodal direction. 

PACS numbers: 71.10.Fd, 71.27.-^a, 74.72.-h, 79.60.-1 



In spite of intense theoretical studies over the two last 
decades, the debate on the microscopic origin of the so- 
called stripe phases characterized by the combined spin 
and charge order is far from closed)^ So far, two main 
scenarios regarding the origin of stripes have been put 
forward. In the weak-coupling scenario, stripe phases 
emerge due to nesting properties of the Fermi surface and 
provide a compromise between the superexchange inter- 
action, which stabilizes the long-range antiferromagnetic 
(AF) order in the parent Mott insulator, and the kinetic 
energy of doped holes. In this scenario spin and charge 
orders occur at the same temperature or charge stripe or- 
der sets in only after spin order has developed. An alter- 
native strong-coupling scenario comes from the Coulomb- 
frustrated phase separation suggesting that stripe forma- 
tion is commonly charge driveui^ In this case the onset of 
charge order appears prior to spin order as the tempera- 
ture is lowered. However, this scenario does not take into 
account spin fluctuations which might prevent the spins 
from ordering at the charge-order temperature Tqo ^ 

Similarly, on the experimental side the origin of 
stripes remained up to now unclear. On the one hand, 
in the most widely studied stripe-ordered compounds 
Lai.6-xNdo.4Sr^Cu04 (LNSCO>i and La2_:,Ba^Cu04^ 
the charge order observed as Bragg maxima at the wave 
vectors q^ = (±47re,0) with e ~ a; for doping x < 1/8 
was found to set in at a higher temperature Tco than 
the spin order signalized below Tgo by the onset of mag- 
netic Bragg peaks at q^ = 7r(l ± 2e, 1), thus supporting 
the strong coupling scenario. On the other hand, the in- 
terval between Tco a-nd Tso was found not to exceed 15 
K. In contrast, in the recent resonant soft x-ray studies 
of the stripe order in Lai.8_2;Euo.2Sr2,Cu04 (LESCO)(^ 
both temperatures were found to be separated by 35 K. 
Those findings clearly exclude the nesting scenario of the 
stripe formation in favor of the electronic origin. 

In this paper we address melting of the metallic stripe 
order and look for a signature of a pure charge-order 
above Tso- We consider the Hubbard model. 



with a hopping amplitude t [t') between the (next-) 
nearest neighbor sites, respectively, and with U = lOi 
standing for the on-site Coulomb interaction and solve 
it within the real-space dynamical mean-field theory 
(RDMFT). It allows one to handle the leading local part 
of dynamic correlations exactly. Although the long-range 
Coulomb interaction might help to stabilize further the 
stripe order, the great success of the previous RDMFT 
studies at T = was a proof that the correct treatment 
of the on-site interaction alone suffices to stabilize ex- 
perimentally observed metallic stripes in LNSCOii Fol- 
lowing Ref. |7|, we decompose the square lattice into L^ 
stripe supercells with jjl = \ . . .Lc orbitals so that the self- 
energy Y,cj{K ,iu!m) and nonintcracting Green's function 
Go{K,iujm) become Lc x Lc matrices, where the wave 
vectors K span the folded Brillouin zone (BZ). Further- 
more, similarly to the standard DMFT,^^ the RDMFT 
approximation neglects the K dependency of the self- 
energy while allowing for its spatial dependency, i.e., 
Yla-iK,iujm) = ^fj.iy,(7ii'-^m)Sfj,i,. Ncxt, the lattice model 
([T]) is mapped onto a set of effective quantum impuri- 
ties subject to a dynamic bath ^o,o-(*'^m) that has to 
be determined sclf-consistcntly. Here we use a Hirsch- 
Fye quantum Monte Carlo (QMC) impurity solver; it not 
only provides access to finite temperatures but also has 
the advantage of simulating an effectively infinite bath. 

In particular, for a given site-diagonal element of 
the bath Green's function Go,ia-{'i'-^m), chosen to im- 
pose an initial stripe configuration, we use Lc times the 
QMC solver to obtain the interacting Green's function 
Qicr{ioJm) and then we extract through the local Dyson 
equation the corresponding self-energy, T^irjiiuim) = 
^oiai^^in) — G^cr (*<^m)- Ncxt, wc compute a new bath 
Green's function using the self-consistency condition in 
a matrix form;^ g'^,^{iujm) = Ga {i^m)Sp,i, + 12a{iujrn), 



with Qa{i^7i 



G^'{K,; 
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all the impurity sites and hence the actual charge density 
at a given site depends effectively on the charge distri- 
bution at other sites in the stripe phase. The proce- 
dure is repeated till convergence is reached. Finally, we 
represent the lattice Green's function, G~ {K, iuirn) = 
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FIG. 1: (Color online) Temperature dependence of the: (a,b) 
local charge density ri;, (c,d) magnetization l^f |, as well as 
(e,f) double occupancies -D; at the inequivalent sites in the 
SC (left) and BC (right) stripe phase at doping x — 1/8. 
Filled symbols at /3f = 50 were obtained with t' — — f/3. 
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i), in the original BZ as follows, 
s*'=("''-"")G.(if,iw™), and by 



rotation to the real-frequency axis via stochastic analyti- 
cal continuation methodfiS extract the A;-resolved single- 
particle spectral function A{k,uj) = —-lni^^ga-{k,uj). 

We begin by showing in Fig. [T]the temperature evolu- 
tion of the local charge ri; = '^^icj^Cj^a) and magnetiza- 
tion \S?\ = ^K'^it ^ "4)1 densities as well as changes in 
double occupancies D^ = (n.^-i-n.^) in the two possible ver- 
tical stripe patterns at doping x = 1/8: (i) site-centered 
(SC) one with enhanced hole density at the nonmagnetic 
domain walls (DW) that separate AF spin domains and 
(ii) bond-centered (BC) one consisting of hole-rich lad- 
ders with a weak ferromagnetic order on the rungs sep- 
arating AF ladders. In fact, using the symmetry of the 
self-energy in the spin sector, it is possible to reduce the 
number of the impurity sites that must be solved down to 
3 (2) sites in the SC (BC) stripe phase, respectively, in a 
magnetic stripe unit cell with Lc ~ 8 atoms. Cuidcd by 
the experimental data one can expect a complex process 
of melting of the stripe order, possibly via an intervening 
phase with a pure charge orderi^ We find instead that 
the spin and charge orders disappear simultaneously at 
/3i ~ 7 and the overall melting is reminiscent of disap- 
pearance of the mean-field stripe phase in the overdoped 
regime^ This follows from Fig. [1] where one observes a 
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FIG. 2: (Color online) Intensity of the: (a) charge C(g) and 
(b) spin S((i) structure factor maxima at the wave vectors 
corresponding to the SC stripe phase with the charge (spin) 
unit cell containing 4 (8) atoms, respectively. For clarity, C(9) 
was multiplied by a factor of 4. 



simultaneous increase in the charge rij and spin S^ or- 
der parameters on decreasing Tj^iH A possible source 
of discrepancy between the numerical and experimental 
data lies in the absence of spatial fluctuations inherent 
to the RDMFT approach. Indeed, the physics of the 
proximity to phase separation in the cupratcs is not cor- 
rectly captured since the later instability corresponds to 
a long-wavelength instability. Alternatively, the melting 
of the stripe order in the layered cuprates might well 
be controlled by the interplane coupling allowing for a 
phase transition to occur at a finite-temperature. In- 
deed, it is not evident that the inter-plane coupling af- 
fects the charge and spin orders in the same way. We 
also note that the enhancement of the double occupan- 
cies in the low T-regime gives insight into the mechanism 
of the stripe formation. The kinetic energy gain acts as a 
driving force behind the stability of charge stripes along 
which doped holes can propagate coherently as seen in 
the corresponding spectral functions (see later). 

To test the robustness of the stripe phase and exclude 
possible stripe meandering into the AF domains we re- 
lease all the constraints and solve self-consistently an 
8x8 stripe supercell at /3t = 8 and 12 where the use 
of the QMC solver albeit time expensive is nevertheless 
feasible. Further, in order to facilitate the analysis, it 
is convenient to compare the intensity of the cusps in 
the charge C{q) = jq ^j. e.'^''^ {no) {ur) and spin S{q) = 

jfJ2r^^'^i^o)i^r) structure factors at the appropriate 
wave vectors q^ = (±7r/2,0) and q^ = (±37r/4, ±7r). 
Again, as shown in Fig. [2] we find that the tempera- 
ture dependence of both structure factors is very similar. 
Therefore, we conclude that relaxing the constraint has 
no overall effect even in the vicinity of the phase transi- 
tion where thermal fiuctuations are important. 

The signature of the stripe order emerges predomi- 
nantly in the low-energy features of the single-particle 
spectral function A{k,Lu). Figs. [3l[a-c) and H] plot this 
quantity as a function of temperature. Let us first con- 
centrate on the nodal direction of the BZ. Here, we note 
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FIG. 3: (Color online) Evolution of the low-energy part of 
the single-particle spectral function A{k, cu) in the SC stripe 
phase: (a) /3t = 12, (b) /3t = 16, and (c) /St = 50; the SC and 
BC stripe phases are dynamically indistinguishable. Panel 
(d) shows A{k,Lo) in the PM phase at /3f = 50. 



that thermally disordered stripes below /3t — 12 (al- 
beit clearly seen in the order parameters in Fig. [Ij leave 
A{k, to) in the nodal direction almost intact. These find- 
ings are consistent with Ref. [ij suggesting that the nodal 
A{k,uj) in the stripe system may arise in the presence 
of disorder or due to fluctuations of stripes facilitating 
penetration of the AF domains by holes. In contrast, a 
further reduction of temperature strengthens stripe or- 
der and modifies its spectral properties. Indeed, already 
at (3t ~ 16 one finds along the nodal direction two sep- 
arated features originating from the single quasiparticle 
(QP) peak. It results in the depletion of A{k,u!) at the 
Fermi level ep- Consequently, as depicted in Fig. [SJa), 
the QP peak seen at ep a.t (3t = 8 in the density of 
states N{uj) is superseded by a pseudogap that opens up 
at /3i = 16 just below sp- Moreover, when moving to- 
wards the 5* = {tt/2, tt/2) point in Fig. SI spectral weight 
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FIG. 4: (Color online) Temperature dependence of the low- 
energy part of the spectral function A{k, lj) in the SC stripe 
phase in the vicinity of: (a-c) the nodal S = {n/2, it/2) and 
(d-f) antinodal X = (tt, 0) points. For comparison, dashed 
line in panels (c) and (f ) shows A{k, lu) in the PM phase. 



is transferred gradually from the lower to the upper QP 
peak whose maximum is far above ep at the S point. 
Hence, in agreement with the previous T = RDMFT 
studies of the stripe order predicting opening a gap for 
charge excitations along the diagonal directiouji it seems 
that at the RDMFT level, well developed metallic stripes 
are incompatible with the observed in La2-xSr2:Cu04 
(LSCO) and LNSCO nodal quasiparticlesJ^ Altogether, 
noting that the measured nodal A{k, w) is more intense in 
LSCO (possible dynamic stripes) than in LNSCO {static 
stripes) and comparing with the temperature dependence 
of A{k, uj) shown in Fig. HI we confirm that the intensity 
of the nodal A(k, co) provides a measure of the stripe dis- 
order. 

Similarly to the case of the nodal direction, it is only 
above /3t = 12 that we observe the formation of a flat 
band near the antinodal X = (tt, 0) point located just be- 
low sp. This flat band can originate from the renormal- 
ization of the effective mass due to the frequency depen- 
dence of the self-energy - as observed in the DMFT stud- 
ies of models of heavy-fermionsi^ - or from its spatial de- 
pendence. To address this issue, we calculate A{k,uj) in 
the uniform paramagnetic (PM) phase at the same dop- 
ing X — 1/8. As depicted in Figs.Etd) andlH neither the 
flat QP band at the X point develops nor the gap around 
the S point opens up in this case. Moreover, the Matsub- 
ara axis mass-renormalization factor i.e., derivative of the 
imaginary part of the self-energy taken at the DW w.r.t. 
the smallest Matsubara frequency — — "g"^ "" U^^ttT 
does not change across the stripe energy scale /3t ~ 16 
[see Fig. EJb)] at which visible changes in A{k,uj) de- 
velop. Hence, it is the site-dependence of the self-energy 
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FIG. 5: (Color online) (a) Temperature dependence of the 
density of states A'^(a;) in the SC stripe phase and (b) imagi- 
nary part of the self-energy at the DW. Dashed line in panel 
(a) shows N{u}) in the PM phase while in panel (b) depicts 
the result of a linear fitting of data points at the smallest 
Matsubara frequency ujm — nT. 



which is responsible for the substantial differences in the 
angle-resolved and integrated spectral functions between 
the PM and stripe phases in the low-T Hmit. 

Finally, let us briefly address the influence of the next- 
nearest-neighbor hopping t' on the melting of the stripe 
order. On the one hand, Fig.[T]shows that the charge and 
spin profiles at our lowest temperature (3t = 50 are only 
weekly affected by a finite t' ~ — 1/3. On the other hand, 
a destabilizing function of t' is seen as a strong reduction 
of Di which signalizes reduced kinetic energy gain. The 
stripe order becomes more fragile and vanishes already 
at f3t ~ 9. While we have again not found a spontaneous 
tendency towards stripe meandering in the 8x8 unit cell, 
one nevertheless observes enhanced discrepancies w.r.t. 
the t' = case indicating the frustrating role of t'. 



In summary, our RDMFT calculations reveal a strong 
temperature dependence of A{k,uj) in the stripe phase. 
The dominant contribution to the electron mass renor- 
malization at the antinodal points originates from the 
spatial dependence of the self-energy. Furthermore, it is 
only well below the transition temperature that a charge 
gap opens in the nodal direction. The melting of the 
stripe phase is mean-field-like since both charge and spin 
modulations vanish simultaneously. Within the RDMFT 
this was shown to be a robust result with respect to 
both t' and the size of the stripe supercell. The dis- 
crepancies with the experimental situation in LESCO, in 
which a two-step melting is observed, can be understood 
in terms of the absence of spatial fluctuations inherent to 
the RDMFT approach. Additional argument comes from 
the recent dynamic cluster approximation studies show- 
ing that the short-range AF correlations in the PM phase 
are sufficient for the occurrence of phase separatioufii 
Hence, it would be interesting to investigate whether the 
tendency towards phase separation can develop into a 
pure charge stripe order seen in LESCO. While other 
terms such as the long-range Coulomb interaction and 
intcrplanc coupling might also be necessary, a systematic 
improvement of the RDMFT aimed at capturing spatial 
fluctuations is presently under progress. 
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